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Recent J = (para) baryonium interpretations of the BES narrow resonance 
data near the e~e + — > pp threshold suggests the existance of ortho baryonium. 
To assist future searches we study J = 1 states, especially vector meson 
leptonic decays, and report RPA calculations for both light and heavy mesons 
using a Coulomb-gauge QCD-inspired model. Since the <^>(1880) is the only 
missing model state, other discovered J = 1 particles in this region are ortho 
baryonium candidates. 



1 Introduction 

Since the November Revolution entailing the discovery of charmonium and uni- 
versal acceptance of quarks, electron-positron collisions have been an effective 
method for novel hadronic production. The narrow widths of the J/tjj and 
subsequently observed T led to their interpretation in terms of new quark 
flavors (see for example [1]). These mesons were quarkonium states which 
predominantly decayed via the electroweak interaction. In a potential quark 
model their leptonic widths were first calculated [2] using 

r ee = r v ^ e - e + = 167ra| M -^2-|-0(O)| 2 (1) 

with q the charge of the quark in electron units, My the resonance mass and 
i/j(0) the wavefunction at the origin. The quantity |^(0)| 2 also appears in the 
matrix element of the hyperfine interaction [1] and therefore its model value for 
the hyperfine splitting, AMi iyp , in meson spectra (pseudoscalar- vector meson 
mass differences) should be simultaneously tested. This also provides a relation 
between the leptonic width and the hyperfine splitting 

1/9 r K , n 4/3 



AMhw= A^^ ] Mv - (2) 
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However, using this result and the known n-p mass difference, the first pre- 
diction [1] for the charmonium hyperfine splitting failed dramatically. Subse- 
quently, quark models, such as Ref. [3] , incorporating a complex, flavor depen- 
dent hyperfine interaction were able to reproduce the hyperfine splittings in 
various circumstances but fine-tuning was still required since these approaches 
did not include chiral symmetry. The key point is that the pion is the Gold- 
stone boson of spontaneously broken chiral symmetry which is responsible for 
its light mass. Indeed more recent and improved chiral approaches [4, 5] find 
that about 70 % of the n-p splitting is due to chiral symmetry. Consequently, 
attributing this large 600 MeV splitting entirely to the hyperfine interaction 
can lead to an inconsistent description of other, less dramatic spectra split- 
tings in baryons, heavy mesons and light meson excited states. 

A more fundamental treatment is provided by field-theoretical quark mod- 
els that dynamically incorporate chiral symmetry. The Nambu and Jona- 
Lasinio model is one example but the signature contact interaction precludes 
describing radial excitations (as well as confinement). A similar problem 
plagues contemporary covariant formulations of the Schwinger-Dyson and 
Bethe-Salpeter equations that generate unphysical states corresponding to 
relative time excitations. A more promising approach, which we have utilized 
[5], is to diagonalize an effective QCD Hamiltonian in the Coulomb-gauge 
using the chiral symmetry preserving Random Phase Approximation [RPA] . 

In this paper we apply the RPA to predict vector meson leptonic widths. 
In addition to further testing our model by confrontation with data we also 
discuss exotic baryonium states. In view of the recent resonance discovered at 
BES [6] with a nucleon-antinucleon bound state interpretation [7, 8] , we com- 
ment on possible vector states in the 1800 MeV region, focusing on the quark 
model's, as yet undiscovered, 0(1880) that naturally fits in our approach. 



2 Vector Meson Self-Energy and Width 

The width, r, of a hadron is inversely related to its lifetime and can be 
obtained from its self-energy, S(p 2 ), by Ira S(M 2 ) = r/2 for a particle with 
mass M and momentum p. For a vector meson the appropriate self-energy 
contribution at the one e~e + loop level can be computed from the diagram 
in Fig. 1. This yields 

where g(p 2 ) is the meson-photon transition form factor and n(p 2 ) is obtained 
from the photon polarization tensor with form governed by gauge invariancc 



n^{p 2 ) = {g^p 2 -p^Pv)n(p 2 ) 



(4) 
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Fig. 1. Feynman diagram for the one loop e~e contribution to the self-energy 
of a meson (solid line) with momeum p. All hadron information is contained in the 
transition form factor g(p). 



Applying standard QED covariant perturbation theory to evaluate the 
photon polarization tensor yields 

n 2 f A Trfa hAh± £>)) (K s 

Note we have neglected the electron's mass since it is a small effect. Related, 
Eq. (5) also represents the one loop contribution in the zero mass limit 

and, as suggested by Table 1, this mass can also be neglected since is 
essentially equal to r ee within experimental error. This conclusion is reinforced 
by being randomly larger or smaller than r ee , since phase space would 
indicate it to be smaller. Further, the decay width for Y — > t~t + is also known 
and is compatible with the widths for [T fi + and e~e + as well. The effect of 
the fi mass should be maximal for the p and uj widths where the momentum of 
the outgoing particles is not as large. The relevant correction (see discussion 
below), Am^/p > , produces a 10 % suppression which is less than experimental 
error. 



Table 1. Observed vector meson masses and leptonic widths. Partial widths esti- 
mated from the central PDG [9] value for the total width. The error quoted, affecting 
the last digit, corresponds to the uncertainty in the branching ratio alone. 

State M(MeV) r ee (keV) r^jkeV) 



p 


770 


6.81(2) 


6.9(4) 




783 


0.60(2) 


0.76(25) 


<t> 


1019 


1.26(2) 


1.22(9) 


J/1l> 


3097 


5.2(1) 


5.1(9) 




3686 


2.1(1) 


2.1(3) 


r 


9460 


1.32(5) 


1.30(3) 


T(2S) 


10023 


0.52(3) 


0.58(9) 



Performing the Dirac traces in Eq. (5) and employing the tensor integral 
relations for a scalar function F 



Felipe J. Llanes-Estrada and Stephen R. Cotanch 

J q,q„F(q,p) A = 9 -f J k 2 A 

+ ^ jm q -p) 2 -p 2 q 2 }Fd*q 
J q,F A = ^ J q-pF d A q 



(6) 
(7) 



one can obtain scalar integral expressions. Since we only require the imaginary 
(finite) part of the U self-energy, we can ignore all real (known to be UV 
divergent) contributions. For example, we suppress the tadpole integral 

f d 4 q <? = r d 4 q 1 

J (2Ky(f+ie)({q+pf+ie) J (2^) 4 (q + pf + ze 

in this calculation (the real, mass shifts from lepton loops are totally negligi- 
ble). After some manipulation 

4e 2 

Im (i n{p 2 )) = — Im (t J) (8) 

where J is the standard massless boson loop integral 

J (27T)4(? 2 +^)((9+p) 2 +^) ' 1 ' 

The qo integral is easily performed in the complex plane and the remaining 
integral can be obtained with the substitution 

2wiS{E - E r ) . 



E-E r -ie 

The result (identical to using Fermi's golden rule) is 

giving the theoretical leptonic width 

r = 2Im W 2 ) = !gV(M 2 ). (10) 

For completeness, the correction to Eq. (10) due to a non-zero lepton mass, 
m, in the polarization tensor is obtained [10] by replacing the explicit fine 
structure constant with 



/ 4m 2 / 2m 2 \ / 4m 4 \ 

a E M — » olem\1 i — ^ 1 1 + J — aEM 1 1 — i ■ [ ' 1 } 
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Fig. 2. Diagrams for the transition form factor g(p): (a) hadron decay and rescat- 
tering before annihilation; (b) direct quark-antiquark annihilation. 



3 RPA Leptonic Widths 

To evaluate r ee in the RPA we need to calculate the transition form factor 
g(p). As indicated in Fig. 2 there are both meson and quark loop contributions. 
In the next section we show that the quark-antiquark annihilation, Fig. 2b, is 
the dominant process with contribution in the RPA given by 

9 = J&^W [^(0)^(0),^] |^) RPA (12) 

where W is the quark field operator and q the (fractional) quark charge in 
electron units. The RPA vector meson creation operator with spin projection 
i is denoted by V 1 ^ and contains creation and annihilation parts 

V*t = X* afl BlDl - Y* p B a D fJ (13) 

with all remaining quantum numbers represented by a and (3. For J PC = 1~~ 
states there are S and D wavefunction components 

™> - TJ^™ 1151 

Xij(k) = ^=-^=1- (k-aP - \<A ia 2 X D {k) (16) 



with color indices a, 6, c and N c = 3. The normalization condition 

reduces to 
1 



(18) 



k 2 dk (\X s (k)\ 2 + \X D (k)\ 2 - |F s (fc)| 2 - \Y D {k)\ 2 ) = 1 . (19) 



The scalar wavefunctions are solutions to coupled integral equations re- 
sulting from the RPA diagonalization of the Coulomb gauge Hamiltonian [5] . 
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A state-of-the-art QCD confining potential with scale (string tension or gluon 
generated mass) connected to /Iqcd was used along with a transverse gluon 
exchange hyperfme interaction, essentially a Yukawa potential, supported by 
recent lattice results [11], with a range parameter chosen to adjust the mass- 
spectra. 

In evaluating Eq. (12) we note that spin is conserved in the vector meson- 
photon transition (diagonal coupling). Rotational invariance permits summing 
over the three spin components if we divide by 3. The resulting RPA g is then 



eq^/MNl f°° k 2 dk 



I 

Jo 



(Xs(k)-Y s (k))(2 + sm4>{k)) 



3 J (2tt)3 
-V2(l-an<f>(k))(X D (k)-Y D (k))\ . (20) 

Note that the RPA leads to an expression of the form given by Eq. (1) 
since r oc g 2 , and in the non-relativistic limit 

g^eqy/^N c J J^*s(k) (21) 

because sin 4>{k), the BCS vacuum angle (solution of the gap equation for the 
same Hamiltonian), approaches unity 

sin m = M ^ ~ 1 (22) 
^ ' ^k 2 + M q (k) 2 

for large running quark mass and the Y function can be neglected in the same 
limit. Finally the Fourier-transformed wavefunction is 

The effective charge factors, q 2 -^>< q 2 >, for the various quark flavors 
necessary for our width calculation are listed in Table 2. This is the origin of 
the order of magnitude suppression of the u> width relative to the p. We work 
in the exact isospin limit. 



Table 2. Quark charge factors (q 2 )- 



State 


Flavor 


(l 2 ) 


r 
p 

LO 


bb 

cc 
ss 

uu-\-dd 

V2 _ 
uu — dd 
V2 


1 

9 

4 
9 

1 

9 

1 

2 

1 

18 
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Combining Eqs. (10) and (20) yields the final expression for the RPA 
widths. The numerical predictions for both meson masses and widths are 
listed in Table 3. The predicted masses are only in rough agreement since 
a fit was not attempted (the quark masses were not fine-tuned). To avoid 
compounding errors, the physical meson masses were used in Eq. (10) except 
for the unknown 0(1880) where the model mass is adopted. Given the large 
errors in the measured lepton widths this should be a minor concern. For 
comparison, two other theoretical calculations [3, 12] are also presented. It is 
noteworthy that our few parameter approach provides a description which is 
comparable to these multi-parameter models. 



Table 3. Experimental and theoretical vector meson widths. Also listed are two 
other calculations, marked (Godfrey and Isgur [3] S-D wave results without mix- 
ing) and 2 (Ebert, Faustov and Galkin [12]). 



State 


Exp. M(MeV) Calculated 


Exp. r e - e +(keV) Calculated Other works 


P 


770 


795 


6.85 


4.7 


6.87 1 


UJ 


782 


795 


0.60(2) 


0.50 


0.85 1 


<t> 


1019 


1005 


1.26(2) 


1.1 


2.75 1 


UJ 


1420 


1420 


0.1(1)* 


0.18 


0.046 1 


p 


1450(25) 


1420 




1.5 


0.45 1 


UJ 


1650(25) 


1620 




0.23 


0.013 1 




1680(25) 


1670 


0.5(10) f 


0.21 


0.27 1 


p 


1700 


1520 


3.5(5)+ 


2.0 


0.14 1 


<t> 




1790 




0.45-0.55 


0.06 1 


J/i> 


3097 


3130 


5.2(1) 


3.8 


9.9 1 , 5.4 2 


1> 


3686 


3681 


2.2(2) 


2.0 


3.3\2.4 2 


i> 


3770 


3695 


0.26(4) 


0.68 


0.10 1 


y> 


4040 


4140 


0.75(15) 


1.7 


1.8 1 


y> 


4160(20) 


4150 


0.77(23) 


0.56 




i> 


4415(6) 


4535 


0.47(10) 


0.26 




r 


9460 


9460 


1.32(5) 


0.45 


1.4\l.3 2 


r 


10023 


9870 


0.52(3) 


0.35 


0.65 1 , 0.5 2 


T 


10365 


9921 


0.5(1)* 


0.0003 


0.45 1 


T 


10580 


10190 


0.32(3) § 


0.29 


0.34 1 



our estimate * fi + fi width (roughly equal) § recently reported by BaBar [13] 



Finally, we can relate our form factor parameter, g, to the standard vector 
meson leptonic decay constant, fy, defined by 

r e - e+ = y«| M My/ y . (24) 
Comparing to Eq. (10) yields 



/ M 3 

9 = efv\Hr- (25) 
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with mass dimension consistent with Eqs. (19) and (20). In representing the 
widths from Ref. [3] we use their fy values and again available physical meson 
masses. Combining the different models it is clear that we have a reasonable 
theoretical understanding of vector meson dileptonic decay. 



4 Mesons with Broad Hadronic Widths 

For vector meson decay to a particle-antiparticle hadron pair, annihilation to 
a virtual photon is possible which also contributes to the total lepton width. 
In this section we will demonstrate that such rescattering corrections are small 
in comparison with the intrinsic quark-antiquark annihilation contribution. A 
typical example of a broad vector meson width is the p — > irir. We assume 
that the true eigenstate of the full QCD Hamiltonian also contains a two-pion 
wavefunction component 

\p) = a\qq) + (3\nn) +■■■ (26) 

with |/3| < 1. We have previously calculated the decay of the qq component 
into e~e + and now proceed to estimate the contribution from the \tttt) term. 

In principle this higher Fock-spacc wavefunction component could be com- 
puted from the same model Hamiltonian employed in the qq two-body sector 
since it is field-theory based. However this is a challenging calculation involv- 
ing four relativistic particles and chiral symmetry. Although a formalism has 
been developed [14] to treat this problem, for an order of magnitude estimate 
we instead implement a simple model wavefunction and adopt the following 
ansatz for the V = p two-pion component 



v l a H P = (M Vl o))\n) 



d 3 k V3fc 4 e abc r wv /2 
(2tt) 3 y/2 M v - 2E k - ir v /2 



2(2tt)3 1 b1 f 
\/2M^Vl k ~ kl ^ ( } 

Here E k = y 'M% + fc 2 , Jy is the total V width and r n7r is the width for 
the 7r + 7r~ decay. The indices a, b and c now denote charge states. This wave- 
function is inspired by the scattering amplitude for h + h~ hadrons coupled to 
angular momentum 1 near a resonance described by the Breit-Wigner ampli- 
tude 

MyThh/2 

aiis) = - s -M v + lMv r v (28) 

or in the V rest frame with E = My + iTy /2 

ai{E) = M v -E- I r v /2 = 2T V -- (29) 
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The normalization is given by 

(n\ \v i a ,v i a 't] \n) 



{Q\Q) 

where the factor N 2 (smaller than unity) is 



§aa §v (2ir) S (0) N (30) 



^ = ^T I dE - ^ ^ (31) 



™ v J 2m ^ {E-M V )*+I*/A 

and represents the probability of finding a pion pair with relative momen- 
tum given by the Breit-Wigner distribution. This determines the coefficient 
above. This normalization also represents the number of pion pairs per p under 
the Breit-Wigner curve. 

The pion pair annihilation to a photon involves the electromagnetic current 

f a (0) = ie ahc 7r bt (0)dV ct (0) (32) 

taken bare as all pion rescattering effects are already included in the Breit- 
Wigner distribution. The pion field at the origin can be expressed in terms of 
the pion momentum operator 



A 0)=/^«< ' . (33, 



d 3 q 

q m ^ q 

The hadronic transition coefficient is then 



where the charge factor is now 1 as both hadrons carry one electron unit and 
spin and isospin are averaged over. The straightforward computation yields 

r k 5 / 2 dk jwg 

9 eq ^Z{2«fM v Jo E 3 k /4 M V - 2E k - zIy/2 ' [ > 

The imaginary part of the integral is finite, but the real part diverges linearly 
with increasing momentum cut-off A. This reflects short range physics and the 
need for a counterterm in the vector current of Eq. (32) . Since the "contact" 
part of the decay has been calculated using the RPA, the counterterm coef- 
ficient (or simply the cutoff in the integral) could be fixed by the difference 
between the physical and RPA lepton widths. But this generates an unnatu- 
rally large cutoff which is not expected for pions with energies comparable to 
the mass V, i.e. pion pairs with energies beyond a GeV from the resonance 
should not appreciably contribute. A third way to regularize the integral is to 
allow for an energy-dependent width that localizes the integrand. In principle 
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Table 4. Contribution to r v _, e - e + from tttv pion rescattering. 



M v {MeV) r v 


(MeV) r nn 


(MeV) 


Af (A — 


M v )/r v 


r ee (keV) 


700 


150 


150 


0.38 


5 


0.068 


700 


150 


150 


0.38 


10 


0.156 


1465 


350 


100 


0.12 


5 


0.028 


1465 


350 


100 


0.12 


10 


0.064 


1700 


240 


80 


0.11 


5 


0.0070 


1700 


240 


80 


0.11 


10 


0.014 



all three ways are equivalent and we therefore use the simple cut-off regular- 
ization with parameter A. The results for the e~e + widths corresponding to 
the three lowest p states are presented in Table 4 for different model parame- 
ters. The momentum cut-off is listed dimensionlessly, (A — My)/7V, in terms 
of the V mass and width. 

Note that the rescattering contribution to the electron-positron widths are 
suppressed relative to direct quark-antiquark annihilation by an order of mag- 
nitude. With this estimate and the present experimental width precision, the 
RPA calculation appears to provide a sufficient description without the need 
for rigorously including hadron final-state interactions. Further, rescattering 
effects for other meson decays can be even smaller. For example, in odd G 
parity w-type meson decays to three or more (odd) particle states, the final 
state mesons are very unlikely to recombine and annihiliatc to a photon. Also 
the 4> an d several other quarkonium states have narrow widths because they 
are near or below hadronic decay thresholds. Finally, the </>(1680) has an uo- 
like decay pattern with dominant K*K decay products which cannot directly 
annihilate into a single photon. 

5 Towards baryonium 

The short-ranged (weak) NN interaction, Vnn, can only produce a single 
bound state, the deuteron. In contrast, for the AW system the stronger V N ft 
interaction obtained from Vnn using G parity can support several bound 
states. However, because of annihilation and also mixing and decay to other 
baryon number zero hadrons (mesons), it has long been believed that bary- 
onium states, if found, would be quite broad. Significantly, BES has re- 
ported [6] a dramatic narrow enhancement in the decay J/ip — > jpp but 
not in J/ip — ► n pp. This indicates the possibility of an I = 0, as opposed to 
1=1, bound state not far below threshold [8]. The quantum numbers of this 
state are likely to be either J PC = h or ++ . If it is a pseudoscalar state, it 
cannot be explained in a qq model since the expected excited n and r\ states 
in this region have already been identified as the 7r(1800) and 7/(1760). 

The baryonium interpretation of this resonance suggests the pseudoscalar 
assignment corresponding to para baryonium (nucleon and antinucleon spins 
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antiparallel) which, in turn, implies a partner, ortho baryonium, with spins 
alligned having total J PC = A similar state is also suggested by the 

observed low-energy proton timelikc form factor behavior as well as by a 
resonance prediction using Vector Meson Dominance [15] just 
below the e~e + — > pp threshold. Further, a state, with mass 1870(10) MeV 
and width 10(5) MeV, has also been reported by the Fenice collaboration 
[16] in e~e + production but awaits confirmation. Another discrepancy with 
expected phase-space is reported [17] in B — > ppK. 



Since baryonium has conventional meson quantum numbers, it is difficult 
to disentangle from qq and 2q-2q states. In particular, quark models predict 
numerous four-quark states that likely overlap forming a hadron continuum. 
The best prospect for observing baryonium states would be if they had nar- 
row widths (possibly due to significant quark recombination) and non-strange 
decay products near the NN threshold. Such states are predicted by the Res- 
onating Group Method [18]. Because the BES and FENICE resonances both 
have widths of order 10 MeV, Coulombic baryonium states bound by just 
QED are unlikely to be observed since these would have very narrow leptonic 
widths given by 



that are of order 10 eV. Instead, we would anticipate one ortho and one para 
baryonium level each with a typical width of order 10 MeV. 

Experimental information above threshold should be regarded more cau- 
tiously. As recently discussed [19], steep peaks in the nucleon timelike form 
factors and J/ip decays could be simple threshold cusp effects. Consequently, 
baryonium claims should focus below threshold on isolated peaks or dips in 
an observable like the unconfirmed dip at 1870 MeV in the total e~e + anni- 
hilation cross section [16]. 



Since ortho baryonium is a vector state it should be readily accessible with 
new high energy and luminosity e~e + colliders, such as the envisioned DA^NE 
upgrade to 2 GeV [20]. In this energy region, the most likely hadron to be 
observed is the excited (f> meson, a combination of 25 and ID waves, with 
mass 1880 MeV predicted by Isgur and Godfrey and 1790 MeV in our RPA 
model. As listed in Table 3, our calculated e~e + width for this state is about 
half a keV. 

We can also phenomenologically estimate this width using known data and 
flavor symmetry. Recall Weinberg's first sum rule 



-a 5 EM M p /(2n) 



(36) 



6 0(1880) 




p* p^e e+ 



e~e 



(37) 
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which appears reasonable [21] since the left and right sides of this equation 
are 1.758 and 1.754 MeV 2 , respectively. Assuming ideal mixing for the u)-(f> 
system, we obtain 

2M p r p ^ e - e+ = 9M r^ e _ e+ = 18M u r u ^ e - e+ (38) 

which is satisfied to within 30 % (10.5, 11.6, 8.4 MeV 2 , respectively). Com- 
bining Eqs. (37) and (38) with the known [9] relative branching ratio 

-^'p(1700)->e-e+-^'p(1700)->2(x7r) 
^p(1700)->all 

we obtain the estimates, -T p (i7oo)^e-e+ = 3.5(0.5) keV, ^i(i880)^e-e+ — 
0.7(3) keV and r , w (i650)->e-e+ — 0.4(2) keV . These are somewhat larger 
but still in agreement with our model calculations. Therefore, in searching 
for J = 1 baryonium candidates below the NN threshold, there may be a 
0(1880) state with r ee about half a keV which should be distinguishable by 
its kaon decay modes. 

7 Outlook 

Recent claims have again raised the question concerning the existence of bary- 
onium states but more solid information is needed near the NN threshold. 
The upgrade of DA<£NE to the 2 GeV range will provide a good opportunity 
for further searches. A possible fruitful strategy would be to focus on e~e+ 
annihilation into an even number of pions in isoscalar 1 states. Resonances, 
although not narrow, are anticipated and there are structure precursors [22] 
from other experiments. By extracting the leptonic widths, the missing <f> state 
below 2 GeV should be found and we predict its lepton width to be around 
half a keV. Its hadron width should be similar to the w(1650), about 200 
MeV. Any other observed vector state is a candidate for ortho baryonium. 

This work is supported by grants FPA 2000-0956, BFM 2002-01003 (MCYT, 
Spain) and DE-FG02-97ER41048 (US Department of Energy) and has been 
commissioned for the "Encuentro de Fisica Fundamental Alberto Galindo" 
meeting at Univ. Complutense, November, 2004- The authors are grateful to 
the workshop organizers and congratulate Professor Galindo on his lubilaeum. 
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